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ABSTRACT 


The two-dimensional orthogonal function expansion technique was used to 
model discrete ocean bathymetric data as a function of cross-range and down- 
range, and to model discrete sound speed-data that is a function of depth and 
down-range. The technique was tested and validated via computer simulation. It 
was shown to have many positive features, but the approach is flawed by the set 
of generating functions that were selected. 

A reflection angle algorithm for arbitrary two-dimensional ocean bottom 
models was tested and validated, and was found to be accurate. Computer 
simulations of acoustic ray interaction with various two-dimensional ocean bottoms 


modeled with the orthogonal function expansion and using the reflection angle 


algorithm were also conducted. 












Accesion For 


NTIS CRA&l 
DTIC TAB 

Unannounced 
Justification 





eee tt Poe er 








By baie coateteaeiet te 


Distribution | 


Avail and] or 
Dist Special 






V1 





L.. INTRODUCTION, oi teice hoo 0k 68 eee tee whee eee 1 
Il. TWO-DIMENSIONAL ORTHOGONAL FUNCTION EXPANSIONS OF 
BATHYMETRIC DATA: dciney cute hed asuiueea eeu te buesacea ees 5 
A. ORTHOGONAL FUNCTION EXPANSION AND THE BASIS 
FUNCTIONS. c¢siscd seat gt aston es acerabeedenuee eee 5 
1. Definition of the Two-Dimensional Orthogonal Function Expansion ..... 5 
2. Derivation of the Basis Functions from the Generating Functions ........ 7 
3. The Generating Functions and the Minimum Mean-Squared Error 
Criterion: 4.2.0:phercascatutds aes 4 boa ea ee eae ae onic 11 
B. COMPUTER SIMULATION RESULTS OF BATHYMETRIC FITS ...... 13 
Ls Exact. Fit: Test Cases o.i5.200 targus tins eis tear settee ee’ 14 
2. Complicated Test Cases. scc.cicastivtcts ee vaeek esp erataer siete 34 
Do SUMMARY (2s saa da eaceeGre yee t ae had ranean eau ele tia ers 51 
lil. THE BATHYMETRIC MODEL APPLIED TO RAY TRACING IN THE RRA 
PALGORTIIM ster oie-bteesee ence tes cease eee aaneae awe eat ae hee 53 
A. REFLECTION ANGLES OFF THE TWO-DIMENSIONAL OCEAN 
BOTTOM MODDED 6465:sc0-tawe Satie weed eotatay ies eeey Es 53 
B. COMPUTER SIMULATION RESULTS .......... 0... cece eee cece eee 58 
IV. TWO-DIMENSIONAL ORTHOGONAL FUNCTION EXPANSIONS OF 
SOUND-SPEED‘DATAS cide oea ce dey tas Oebsle eo dae eee oh 75 
A. ORTHOGONAL FUNCTION EXPANSION AND THE BASIS 
FUNCTIONS 3c 4.cseacdt eae ee ne eee ee he ee ee eke ease eens 75 
1. Definition of the Two-Dimensional Orthogonal Function Expansion ..... 1] 
2. The Basis Functions from the Generating Functions .................4. 77 
B. COMPUTER SIMULATION RESULTS OF SOUND SPEED FITS ....... 79 
V. CONCLUSIONS AND RECOMMENDATIONS .......... 00s eee eee eee, 97 
LIST OF REFERENCES: cca. neeestvensia eee io atanaves Sdeec eminent 101 
INITIAL DISTRIBUTION: EIS E ‘os ccnriadietereeries seater kneatua oss scceser 103 





TABLE OF CONTENTS 


vil 





I. INTRODUCTION 


In the area of underwater acoustics, the development of an accurate, fast, and 
simple computer model to calculate ray paths and sound-pressure levels is an important 
first step for modeling and studying sound propagation in the ocean. The RRA (recursive 
ray acoustics) algorithm developed by Professor L. Ziomek [Ref.1,2] provides these 
features. It requires only discrete input data of the bottom depth and the sound speed for 
the ocean being modeled, and uses orthogonal function expansions to fit the data in a 
minimum mean-squared error (MMSE) sense using orthogonal polynomials. However, 
until now, the concept of orthogonal function expansions had only been implemented and 
tested for one-dimensional data in the RRA algorithm, that is, bottom depth as a function 
of down-range and the speed of sound as a function of depth. 

The intent of this thesis is to show that two-dimensional orthogonal function 
expansions can be used to fit discrete data as a function of two independent variables, that 
is, bottom depth as a function of cross-range and down-range, and the speed of sound as a 
function of depth and down-range. This technique is implemented in the RRA algorithm 
to provide more realistic models of the ocean bottom and the speed of sound. 

The direction of propagation of a sound ray reflected off the ocean bottom is 
calculated based on the unit normal vector to the ocean bottom at the point of incidence. 
The unit normal vector is calculated based on the derivatives of the ocean bottom with 
respect to down-range and cross-range. 

For a one-dimensional model of ocean bottom depth (as a function of down- 
range) the angle of reflection is computed only from the derivative of bottom depth with 
respect to down-range. The one-dimensional model provides no cross-range information; 
the derivative of bottom depth with respect to cross-range is zero. A two-dimensional 
model of ocean bottom depth (as a function of cross-range and down-range) would 
provide derivatives of the ocean bottom with respect to both cross-range and down-range. 
A two-dimensional math model will provide components in the down-range and the 


cross-range direction, and the reflected ray will travel down-range and cross-range. By 








modeling the ocean bottom with a two-dimensional orthogonal function expansion, an 
accurate model of the ocean bottom will be incorporated into the RRA algorithm. 

The sound-speed model currently implemented uses only one-dimensional data to 
produce a sound-speed vs. depth profile that is used in a zone (a particular area of the 
ocean). A two-dimensional orthogonal function expansion results in a model of the 
sound-speed as a function of depth and down-range that could be used to fit multiple 
sound-speed profiles (SSP) that are taken at different down-range locations. 

Chapter II discusses the two-dimensional orthogonal function expansion used to 
fit bathymetric data as a function of down-range and cross-range. The generating 
function that was chosen, and the algorithm used to orthogonalize the set of generating 
functions using the Gram-Schmidt procedure is detailed. Computer simulation results are 
shown that demonstrate the accuracy and the applicability of this method. The result is a 
smooth surface fit of the ocean bottom obtained from discrete data. 

A shortcoming of the generating functions chosen is that cross-range and down- 
range values have an equal effect on the estimated value of the ocean bottom depth. In 
some cases, ocean bottom depth may have little or no dependence on one of the 
parameters. Independent weighting of the parameters will correct this shortcoming for all 
planar ocean bottom cases. This thesis will demonstrate the need to weight parameters 
independently in these cases, and show some examples where weighting is applied. For 
more complicated ocean bottom geometries, a simple weighting scheme will not correct 
the problem created by the choice of generating functions. Alternative choices for the set 
of generating functions used in the orthogonal function expansion is discussed in Chapter 
V, and is offered as an area for future research. 

Chapter III presents the implementation of the ocean bottom models developed in 
Chapter II into the RRA Algorithm FORTRAN ray tracing program. The orthogonal 
function expansion technique allows accurate calculations of ray reflections off the ocean 
bottom. The algorithm used to calculate accurate reflection angles for sound rays will be 


presented. Using the unit normal vector to the ocean bottom at the point of incidence, the 








direction of propagation of a reflected ray is obtained using vector calculus. Computer 
simulation results showing incident and reflected rays propagated with the RRA 
algorithm are presented. 

Chapter IV discusses the application of the two-dimensional orthogonal function 
expansion technique to fit sound-speed data as a function of both depth and down-range. 
The procedure is the same to fit both bathymetric and sound-speed data. Computer 
simulation results are shown that demonstrate the smooth fit of the discrete data. The 
results show that one math model can fit more than one sound-speed profile (SSP) taken 
at different down range locations. 

Chapter V discusses the problems encountered with the two-dimensional 
orthogonal function expansion technique used to model bathymetric and sound-speed 


data, states conclusions as well as recommendations for future research. 











Il. TWO-DIMENSIONAL ORTHOGONAL FUNCTION 
EXPANSIONS OF BATHYMETRIC DATA 


A. ORTHOGONAL FUNCTION EXPANSION AND THE BASIS FUNCTIONS 


Ocean bathymetry data is commonly measured and stored as a set of discrete 
measurements of ocean bottom depth at distinct, sampled locations. Matrices of data 
have been generated by the U.S. Navy to chart ocean bottom depth, and provide bottom 
contour navigation. This type of data is ideal for developing a model of the ocean bottom 
that can be used in the RRA (Recursive Ray Acoustics) algorithm. In ray acoustics, the 
model of the ocean bottom is used to calculate ray reflections at the point of incidence to 
the ocean bottom. The model of the ocean bottom must provide a smooth fit to the 
_ discrete data, and it must be capable of providing first and second order derivatives of the 
ocean bottom at the point of incidence. The orthogonal function expansion technique 
provides these features, and it has been demonstrated as an effective method when ocean 
bottom depth is a function of one variable, that 1s, down-range [Ref. 2]. 

1. Definition of the Two-Dimensional Orthogonal Function Expansion 

The two-dimensional orthogonal! function expansion used to fit a smooth surface 
toa rectangular matrix of ocean bathymetric data as a function of cross-range and down- 


range is given by 


3,(x2) = 4,0, (x2), (2.1) 


n=0 


where the coefficient 


b, = (ys, (x; fy Ps, (x, Sy )) > 


M,, Ms, 


= > D> Vo, (Xp 2%): (Xp) n=0,1,--°,N, , (2.2) 


i=] j=l 





represents the inner product of the measured discrete ocean bathymetric data y, (x, ,2, ), 
for i=1,2,...,M, and j= L2eceM; , with the set of orthonormal basis functions 
QP, (x,z), N=0,L...,N, » evaluated at the discrete cross-range (x, ) and down-range (z, ) 
values corresponding to the discrete ocean bathymetric data. The order of the expansion 
is represented by N, . The discrete data must be entered in a matrix, where M, is the 


total number of cross-range rows and M, is the total number of down-range columns. 


It is important to note that although the data is entered in a matrix, the matrix is, in 
general, rectangular, not square, and the data is, in general, unevenly spaced in both the x 
and z directions. Figure 1 shows an example of unevenly spaced data entered in a 


rectangular matrix. Each data point of the matrix (i,/) has a cross-range (x, ), down-range 


(z,), and depth y, (x,,z,) value associated with it. 
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Figure 1. Discrete bathymetric data, unevenly spaced, entered in a rectangular matrix. 


2. Derivation of the Basis Functions from the Generating Functions 

The set of orthonormal basis functions @ b, (x,z), n=0,1,...,N,, were determined 
by using the Gram-Schmidt procedure [Ref. 3] on the set of generating functions 
f,(x,z) =(x+z)", n=0,1,...,N,. The choice of this particular set of generating 
functions will be discussed later. The general equation of the Gram-Schmidt 


Orthogonalization Procedure is given by 


(F, (252), 9», (%2)) 
(p,, (x,Z), Ps, (x,z)) 
(F,(%52),Ps,, (*2)) 


oo, (2,2), 2 = 0,1,2,.... (2.3) 
(,,, (x,Z), Ps, (x,z)) 


P., (%2) = f,(%2)— 0, (2)-. 


The Gram-Schmidt procedure to find 9, ,9,,and g, is shown in Equations (2.4) to 


(2.16). The development of the higher order basis functions, for n = 3,4,...,N,, is not 
shown since the procedure to expand Equation (2.3) is similar to that demonstrated with 
the first three basis functions. 

The Oth-order basis function is given by 


PD, (X52) = fy(x,z) = (x + z) =1. (2.4) 
The /st-order basis function is given by 


(Ff, (2,2), (%,2)) 


s2)s (2.5 
(0, (x,2)., (x,2)) Pr, (x,Z) ) 


P ;, (x,2) oe f (x,Z) = 








where the inner products on the right-hand side of Equation (2.5) are given by 


M,. Mz, 


(Ff, (242).0s, (%2)) = > Ay .2y 0%, Hp 2y) > 


i=l j=l 


My, My, 


=D dt, +2)", 


i=l j=l 


My. Mp, 


ES +4, 2.6) 


i=l j=l 
and 


M, by My, 


Eo, =e (,, (x,Z), 9, (x,z)) = » YO, (x, Zi ), (x, »Z;) : 


i=) j=l 


M,, M,, 


=>'D1-1=M, M, , (2.7) 


i=l j=l 


is the energy contained in g, (x,z). Note that the ratio of inner products that appears on 


the right-hand side of Equation (2.5) will show up in later equations. It needs to be 


evaluated during the computation of the /st-order basis function and is saved as 


7, = iGO 2) _(AG2)Ps 2) 


= = 1 —____—___—_ , (2.8) 
(py, (x, Z), P», (x, z)) Eo, | 
Finally, the /st-order basis function becomes 
D,(x,Z)=x+Z—-Wy . (2.9) 





Similarly, the 2nd-order basis function is given by [see Equation (2.3)] 


(f, (X52), 9, (x,z)] 
(Ds, (X,Z),D,, (x,z)) 


_ (F220, 052)) 
(9s, (*,2).Pz, (%,2)) 


P>, (x,Z) = I, (x,2Z) - P », (x,2) 


g,(%2), (2.10) 


where the inner products in the numerators on the right-hand side of Equation (2.10) are 


given by 


My, Ms, 


(F, (52). (:2)) = YD. fy (%p 20%, (%) 2p), 


i=} j=) 


Ms, M,, 


=> > (xy +z)’ . (2.11) 


i=l j=l 


and 


Ms, Ms, 


(f, (x,2), Pr, (x,z)) = py vf, (x; Zi Ms, (x, <i ) 9 


i=] j=l 


My, M, My, Ms, 


= LLG, +2) - Woo D(a +24) (2.12) 


ix] j=l i=} j=l 


The energy contained ing, (x,z)is given by Equation (2.7) and the energy contained — 


in @, (x,Z) is given by 


My, Ms, 


Ey, = (Pp, (%22)Px, (%52)) = >, 2, Pa, X92), %y02y) » 


i=1 j=l 


M,. My, 


=F Son Oy .5)F. 2.13) 


i=] j=l 


Equation (2.13) is evaluated by using Equation (2.9) for the /st-order basis function, 
9,,(x,z). The ratios of inner products that appear on the right-hand side of Equation 


(2.10) are saved for use in later equations as 


(5, 2,0, (2) E,, -_— 
and 
(F,(%,2),01,(%2))  (f:(%,2). 04, (%2)) 
"To Gag) Ey (2.15) 
Finally, the 2nd-order basis function is given by 
Py, (X52) = (x +z)’ —W,, —W,,[x+z-W] . (2.16) 


As mentioned earlier, the equations for the higher-order basis functions, for 
n=3,4,...,N,, are found by expanding Equation (2.3) in the same manner. It is 
important to note from Equations (2.4) to (2.16) that all the input needed to determine the 
basis functions comes from the discrete data points of the ocean bathymetry. This 








technique is ideally suited for implementation by a computer program. Additionally, the 
order of the orthogonal function expansion, N, , can be increased by simply computing 


the next higher order coefficient 5, and orthonormal basis function ¢, ; all previous 


coefficients and basis functions do not have to be recomputed. In fact, many of the 
computations needed for the lower order coefficients and basis functions can be saved 
and used in the computation of the higher order coefficients and basis functions. The 


choice of N, will be discussed in the next section concerning the minimum mean- 


squared error (MMSE) criterion for the best fit of the data. 

Finally, the data for this technique do not have to be evenly spaced. The ability to 
use unevenly spaced data allows for data sets with slight perturbations in the down-range 
or cross-range values. Real data will contain randomness. For example, set and drift of a 
ship collecting bathymetric data could introduce such perturbations in down-range and 
cross-range of measured data, but there would be no effect on the estimated value of 
ocean depth. 


3. The Generating Functions and the Minimum Mean-Squared Error 
Criterion 
Theoretically, any finite set of functions may be transformed into a set of 


orthonormal! basis functions by using the Gram-Schmidt procedure, and a MMSE fit to a 
finite set of data can be obtained using Equation (2.1) [Ref. 3]. Basis functions that are 
based on sinusoidal functions are frequently used in signal processing to model natural 
occurring phenomena with a finite Fourier series expansion. A previous effort to model 
the ocean bottom using a finite Spatial Fourier Series (SFS) expansion was attempted in 
1991 [Ref. 4]. The SFS technique encountered difficulty due to the oscillations that are 
inherent in using a finite set of sinusoids to fit the ocean bottom. Additionally, the 
estimation of the first and second order derivatives at the point of incidence of a ray 
reflecting off the ocean bottom proved to be inaccurate due to the oscillatory fit. The use 


of a set of generating functions with increasing powers of the independent variable, 








f(z)=z", n=0,1,...,M,, was used successfully for one-dimensional orthogonal 


function expansions of ocean bathymetric data [Ref. 2]. 

The choice of using increasing powers of the sum of the cross-range and down- 
range variables as the set of generating functions was used in this thesis research, that is, 
f,(x,z) =(x+z)", n=0,1...,N,. The lower order functions, for n = 0,1, and 2, fit the 
most common trends in the ocean bathymetry without introducing oscillations. For 
n = 0, the orthogonal function expansion fits a constant depth ocean bottom. For n=1, 
the orthogonal function expansion fits an up-sloping or down-sloping ocean bottom in 
both the x and z directions. For »=2, the orthogonal function expansion fits a slope 
with increased curvature. The higher orders tend to fit ocean bottoms with many and 
various contour changes. In computer simulation, the highest order fit used is the 
designer’s choice. Section B of Chapter II will present results with the highest order fit 
being 7th-order, that is, N,=7. The choice of the highest order fit is a trade-off 
between exactly fitting data (interpolation) versus providing a smooth fit to the trend of 
the data. Higher order, or a higher value of N, , will fit contour changes more exactly. 
That is, the estimated fit will more closely match the input data. A lower order fits the 
discrete data by smoothing through contour changes in the bathymetric data that could be 
due to a rough bottom or measurement error. Data is corrupted by measurement noise 
and by a lack of precision in the measurement equipment. Limiting the highest order to 
N, =7 provides a smooth fit when the number of data points is greater than eight points. 

If the discrete data is fit as accurately with N, <7 as with N, =7, then the use 
of the higher order terms is unnecessary. Lower order fits are used based on a minimum > 
mean-squared error criterion [Ref. 2,3]. The order of the fit is determined by computing 


the mean-squared error ( MSE ) for each order fit, NV, = 0,1,...,7, as follows: 


MSE, =E, -E, . (2.17) 


12 





where 


M,, Ms, 


Ey, = (Vs, (%y 2%, (% 2p) = 2a, Fy s2y)? (2.18) 


i=) j=l 


is the energy contained in the measured bathymetric data, and 


My. My, 


Es, = (J, (x; 96 No (Xj 92, )) > ey, (x, oro (2.19) 


i=l j=) 


is the energy contained in the bathymetric estimate. The value of N, that produces the 


lowest MSE becomes the order of expansion used to reconstruct the bathymetric data. 
B. COMPUTER SIMULATION RESULTS OF BATHYMETRIC FITS 


The technique of orthogonal function expansion to model ocean bathymetry that 
was derived in Section A of Chapter II was implemented in the RRA Algorithm 
FORTRAN ray tracing program. The results of various test cases are presented 
graphically to demonstrate both the successes and the shortcomings of this method of 
ocean bathymetry modeling. The graphs, presented in Figures 2 to 23, show the discrete 
bathymetric data points as large black dots (e). The éstimated fit to the discrete data 
computed by the orthogonal function expansion is presented as a mesh of lines to 
represent the two-dimensional, continuous surface. As can be seen from these figures, 
the surfaces, in general, do not pass through all of the data points. Instead, the curves fit 
the trend of the data. Computer simulation of this technique did reveal shortcomings. 
The test cases presented will show examples of the problems encountered, discuss the 
reasons for the problems, and suggest areas of future research to explore methods to 


correct these shortcomings. 
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Each plot contains a legend with the following information: 
(1) ‘CASE:’, a description of the test case 
(2) ‘SIGMAYB:’, the value of the standard deviation, (co, ), of a randomly 


generated zero-mean gaussian noise added to each data point 


(3) the order of the fit, (NV, ), based on the minimum mean-squared error 


criterion 
(4) ‘MSEO:’ to ‘MSE7:’, the mean-squared error for N, = 0,1,2,...,7, computed 
using Equation (2.17) 

The test cases are grouped in two categories. Oth-order and Ist-order fits are 
presented to show that exact fits are obtained in simple cases. The ability to handle 
unevenly spaced data is presented. Also, a problem with the choice of the generating 
function is discussed. Then, ocean bottoms with a more complicated or arbitrary shape 
are presented to demonstrate fitting the trend of the data. The shortcomings of the 
generating function that was chosen is discussed further. 

1. Exact Fit Test Cases 

Figures 2 and 3 demonstrate the fit of the most basic test case, a flat ocean bottom 
with a constant depth. The exact equation used to generate the 8x10 rectangular 
matrices of data for Figures 2 and 3 is given by | 


y, (x,z) = 1300. (2.20) 


Figure 2 shows the discrete data points entered at evenly spaced cross-range and 
down-range locations. Figure 3 shows that the data may be entered at unevenly spaced 
coordinates without affecting the fit. In both cases the fit is an exact Oth-order fit with the 
orthogonal function expansion returning a value of 1300 meters for any combination of 
cross range and down range. The Ork-order mean-squared error, MSEO, is 9.09 x 107’ 


and not exactly equal to zero due to round-off error in the computer simulation which 
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CASE: 8X10 FLAT OCEAN BOTTOM 
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Figure 2. Flat ocean bottom with evenly spaced data. 
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uses double precision numbers. The values of the higher order mean-squared errors, 
MSE1 to MSE7, are equal to MSEO. This shows that no additional accuracy is obtained 
for NV, >0. 


Figures 4, 5, and 6 demonstrate exact /st-order fits. The exact equation used to 


generate the 7 x 8 rectangular matrices of data for Figures 4 and 5 is given by 
l 
,Z) = 1500 -——(x +z), 2.21 
ys (x52) ico z) (2.21) 


which represents a planar ocean bottom that slopes upward at the same rate with respect 
to cross-range and down-range. The exact equation used to generate the 8x8 


rectangular matrix of data for Figure 6 is given by 
y,(x,2Z) = 1350+ (x +2), (2.22) 


which represents a planar ocean bottom that slopes downward at the same rate with 
respect to cross-range and down-range. In Figures 4 and 6, the orthogonal function 


expansion fits the data with an exact first order fit. For N, =1, the orthogonal function 


expansion is given by 
Vp (X,Z) = by +B, (x+2—Wy) = (by — 5H.) + 5, (x +2), (2.23) 


which exactly matches the form of Equations (2.21) and (2.22). Note that the choice of 
generating functions used for this technique, f,(x,z)=(x+z)”", for n=0,1,...,7, will 
result in an orthogonal function expansion based upon the sum of x and z. Since x and z 


are multiplied by the constant 5,, the same slope with respect to cross range and down 


range will result. 
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CASE: 7X8 UP-SLOPE IN BOTH X AND Z DIRECTIONS 

SIGMAYB: 0.0M  1ST-ORDER MMSE SURFACE FIT 

MSEO: 1.295x10° MSE2: 2.645*10” MSE4: 2.645x10° MSE6: 2.645%10° 
MSE1L: 2.645x10° MSE3: 2.645%10° MSES: 2.645x10° MSE7: 2.645*10" 


Figure 4. Planar ocean bottom sloping upward in both the 
x and z directions with evenly spaced data. 
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CASE: 7X8 UP-SLOPE IN BOTH X AND Z, UNEVENLY SPACED DATA 

SIGMAYB: O.OM 7TH-ORDER MMSE SURFACE FIT 

MSEQ: 1.237%10° MSE2: 4.415*10° MSE4: 4.334™10~ MSE6: 4.29410" 
MSE1: 4.805107 MSE3: 4.341%107 MSES: 4.303%10° MSE7: 4,238x10" 


Figure 5. Planar ocean bottom sloping upward in both the 
x and z directions with unevenly spaced data. 
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Figure 5 demonstrates that the data may be unevenly spaced. Note that the fit 
indicated in Figure 5 is 7th-order and not /st-order. Equation (2.21) was used to generate 
the discrete data entered for the computer simulation shown in Figure 5, however, only 
five significant digits were entered. Round-off error results from having only five 
significant digits in the discrete bathymetric data. For the orthogonal function expansion, 
the double-precision computations that evaluate the best fit of the data based on the 
mean-squared error criterion found that the 7th-order fit was slightly more accurate. The 
st-order mean-squared error is of the same order of magnitude as the 7th-order mean- 


squared error. The orthogonal function expansion coefficients, 


b, = 8.82963 x 10° , b, = -1.11226 x10’, 
b, = 6.24349 x 10° , b, = —2.73240x107 , 
b, = 8.38137x10™ , b, = -1.74958x 107 , 
b, = -9.30174 x 10° , b, = 2.38320x10" , 


show that the fit is essentially /st-order. The Oth- and Jst-order coefficients, 5, and 6,, 
are five orders of magnitude larger than the higher order coefficients, 5, ,5,,... ,b, . In this 
test case, the higher order terms of the orthogonal function expansion will affect only the 
fourth or fifth significant digit of the ocean bathymetry evaluated using the orthogonal 
function expansion. 

Figures 7, 8, and 9 show the results of fitting discrete bathymetric data that shoals, 
or has an upward slope, in only the down-range (z) direction. The data is constant with 
respect to variations in the cross-range (x) direction. The exact equation used to generate 


the discrete data for Figures 7, 8, and 9 is given by 


y,(x,z) = 550- a2. (2.24) 
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Figure 7 shows the results of the orthogonal function expansion as detailed in 
Section A of Chapter II. The fit that is presented in Figure 7 is the best fit to the data 
when using the generating function given in Section A of Chapter II. The corner of the 
surface fit at the maximum cross-range and minimum down-range location (x = 1.250km, 
z=0.0km) is above the discrete bathymetric data points; the corner of the surface fit at 
the minimum cross-range and maximum down-range location (x = 0.0km, z= 1.050 km) 
is below the discrete bathymetric data points. The fitted surface approximates a /st-order 
planar ocean bottom model that slopes upward equally in both the x and z directions. 
Therefore, the fitted plane is above half of the discrete bathymetric data points in the 
lower half of Figure 7, and it is below half of the discrete bathymetric data points in the 
upper half of Figure 7. The best fit is limited in its accuracy due to the fact that x and z 
are summed in every term of the orthogonal function expansion. The orthogonal function 
expansion does not treat the x and z dependence separately, and this is the major 
shortcoming of this technique. 

In order to correct this shortcoming, an independent weighting with respect to x 
and an independent weighting with respect to z need to be added to the orthogonal 
function expansion. With the addition of weighting factors to Equation (2.1), the 


orthogonal function expansion becomes 


Ny 
Y, (xz) = > 5.9%, (W,x,W,z) , (2.25) 


n=0 


where W, and W, are weighting factors to introduce independence in x and in z. For the 
exact /st-order cases presented in Figures 7 through 15, the weighting factors are 
constants in order to independently match the rate of change of the ocean bottom depth 
with respect to x and to z. In the following section, for more complicated ocean bottom 


geometries, the weighting factors would need to be functions of the independent 
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CASE: 6X8 UP-SLOPE IN Z DIRECTION CNLY, NO WEIGHTING 


SIGMAYB: O.O0M 7TH-ORDER MMSE SURFACE FIT 


MSEO: 1.575x10° MSE2: 9.557«10' MSE4: 9.452%10° MSE6: 9.41910" 
MSE1: 9.557x10* MSE3: 9.452«10° MSES: 9.419x10° MSE7: 9.384%10° 


Figure 7. Planar ocean bottom sloping upward in the z direction only with no variation 
in the x direction, modeled without weighting x and z separately, with evenly spaced data. 
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variables. The development of an algorithm to determine the weighting factors, W, and 
W., from the discrete data is not presented as part of this thesis and is left as an area of 
future research. The choices of W, and W, in Figures 7 through 15 are based on a priori 
knowledge of the signs (positive or negative) of the slopes of the data with respect to 
cross-range (x) and down-range (2). | 

Figure 8 presents the same test case presented in Figure 7, except that weighting 
factors, W, and W., were added as shown in Equation (2.25) with W, =0 and W, =1. 
The result is an exact /st-order fit. Figure 9 presents the same test case with the discrete 
bathymetric data entered at unevenly spaced locations. As discussed concerning Figure 5, 
the fit determined by the orthogonal function expansion is 7th-order due to using only 
five significant digits for the discrete bathymetric data. The overall fit is /st-order; the 
higher order terms of the orthogonal function expansion only affect the lower significant 
digits of the estimated value of the ocean bottom. 

Figures 10 and 11 show the fit to discrete bathymetric data that has a downward 
slope in only the down-range (z) direction. The data is constant with respect to variations 
in the cross-range (x) direction. The exact equation used to generate the discrete data for 


Figures 10 and 11 is given by 
! 
Y, (x,Z) = 1700+ 70° ‘ (2.26) 


Figure 10 shows the result of the orthogonal function expansion as detailed in Section A 
of Chapter II. Figure 11 presents the same test case presented in Figure 10, except that 
weighting factors, W, and W,, were added as.shown in Equation (2.25) with W, =0 and 
W=1. 

Figures 12 and 13 show the results of fitting discrete bathymetric data that has 
opposite slopes with respect to cross-range and down-range, that is, an upward slope in 


the down-range (z) direction and a downward slope in the cross-range (x) direction. The 
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CASE: 6X8 UP-SLOPE IN Z DIRECTION ONLY, WITH WEIGHTING 

SIGMAYB: 0.0M  1ST-ORDER MMSE SURFACE FIT 

MSEO: 1.575x10° MSE2: 1.094%10° MSE4: 1.094%10° MSE6: 1.094x10™ 
MSE1: 1.094x10° MSE3: 1.094x10° MSES: 1.094%10° MSE7: 1.094x10™ 


Figure 8. Planar ocean bottom sloping upward in the z direction only with no variation 
in the x direction, using W, = 0 and W. =1, with evenly spaced data. _ 
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CASE: 6X8 UP-SLOPE IN Z ONLY, WITH UNEVENLY SPACED DATA 

SIGMAYB: O0.0M 7TH-ORDER MMSE SURFACE FIT 

MSFO: 1.523x10° MSE2: 3.621x107 MSE4: 3.522«10° MSE6: 2.803*10- 
MSEL: 3.740x10% MSE3: 3.523x10° MSES: 3.149%10° MSE7: 2.803x10" 


Figure 9. Planar ocean bottom sloping upward in the z direction only with no variation 
in the x direction, using W, = 0 and W, = 1, with unevenly spaced data. 
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CASE: 8X8 DOWN-SLOPE IN Z DIRECTION ONLY, NO WEIGHTING 


SIGMAYB: O0O.QM : 1ST-ORDER MMSE SURFACE FIT 


Figure 10. Planar ocean bottom sloping downward in the z direction only with no 
variation in the x direction, modeled without weighting x and z separately, with evenly 
aced data. 
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Figure 11. Planar ocean bottom sloping downward in the z direction only with no 
ariation in the x direction, using W, = 0 and W, =1, with evenly spaced data. 
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slopes are equal in magnitude, but opposite in sign. The exact equation used to generate 


Figure 12 shows the result of using the orthogonal function expansion as detailed in 
Section A of Chapter II. Figure 13 presents the same test case presented in Figure 12, 
except that weighting factors, W, and W,, were added as shown in Equation (2.25) with 
W. =—-1 and W, =1. 

Figures 14 and 15 show the results of fitting a rectangular matrix of 
discrete bathymetric data that has opposite slopes with respect to cross-range and down- 
range, that is, a downward slope in the down-range (z) direction and an upward slope in 
the cross-range (x) direction. The slopes are equal in magnitude, but opposite in sign. 
The exact equation used to generate the discrete bathymetric data for Figures 14 and 15 is 
given by 


y,(x,Z) = 1300-— (x —Z). (2.28) 


Figure 14 shows the result of using the orthogonal function expansion as detailed in 
Section A of Chapter II. Figure 15 presents the same test case presented in Figure 14, 


except that weighting factors, W, and W,, were added as shown in Equation (2.25) with 


W. =-1 and W, =1. 

Figures 2 through 15 have shown that the method of orthogonal function 
expansion accurately models discrete bathymetric data in the most basic test cases of 
planar ocean bottoms. From discrete bathymetric data, a continuous model is produced. 


The discrete bathymetric data does not need to be evenly spaced. Errors in the discrete 
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CASE: 7X8 DOWN-SLOPE IN X, UP-SLOPE IN Z, NO WEIGHTING 

SIGMAYB: O.0M 7TH-ORDER MMSE SURFACE FIT 

MSEO: 1.295«10° MSE2: 1.271»%10° MSE4: 1.267%10° MSE6: 1.265*10° 
MSE1: 1.271%10° MSE3: 1.267%10° MSES: 1.265*10° MSE7: 1.263x10° 


Figure 12. Planar ocean bottom sloping upward in the z direction and downward in the 
x direction, modeled without weighting x and z separately, with evenly spaced data. 
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CASE: 7X8 DOKN-SLOPE IN X, UP-SLOPE IN Z, WITH WEIGHTING 
SIGMAYB: O.GM  1ST-ORDER MMSE SURFACE FIT 

MSEQ: 1.295x10° MSE2: 3.353x10° MSE4: 3.353x10” MSE6: 3.353™107 
MSEL: 3.35310" MSE3: 3.353x10° MSES: 3,353x10° MSE7: 3.353%10" 


Figure 13. Planar ocean bottom sloping upward in the z direction and downward inthe 


x direction, using W, = —1 and W, = 1, with evenly spaced data. 
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CASE: 5X6 DOWN-SLOPE IN Z, UP-SLOPE IN X, NO WEIGHTING 

SIGMAYB: 0.0M 7TH-ORDER MMSE SURFACE FIT 

MSEO: 1.092%10° MSE2: 1.082%10° MSE4: 1.078%10° MSEG: 1.07710" 
MSE1: 1.082»%10° MSE3: 1.078x10° MSES: 1.077%10° MSE7: 1.075%10° 


Figure 14. Planar ocean bottom sloping downward in the z 
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x direction, modeled without weighting x and z separately, 
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Figure 15. Planar ocean bottom sloping do 
x direction, using W, 
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bathymetric data, such as round-off error, do not affect the dominant results of a /st-order 
fit. The only major shortcoming is the need of an algorithm to weight x and z 
independently. 

2. Complicated Test Cases 

Figures 16(a), 16(b), and 17 show test cases that were presented in the previous 
section, except that now the discrete ocean bathymetric data has been corrupted by a zero- 


mean random gaussian noise with standard deviation, o, , that has been added to each 
point. Note that the plots refer to the value ofa, as ‘SIGMAYB:’ in the legend of the 


plot. Only the value of the depth of each point has been corrupted. The cross-range and 
down-range values have been left at evenly spaced coordinates so that the (x,z) coordinate 
of each discrete data point corresponds to the intersection of the lines in the mesh of the 
estimated ocean bottom. 

Figures 16(a) and 16(b) show the flat ocean bottom that was presented in Figure 
2. Equation (2.20) has been modified as follows: 


y,(%,z) = 1300+ 0, N(0,1), (2.29) 


where N(0,1) represents a zero-mean gaussian (normal) random variable with a unit 


standard deviation obtained from a gaussian random number generator. For Figure 16(a), 


o,, is 13 meters, a standard deviation that is 1% of the value of the uncorrupted ocean 


bottom data. For Figure 16(b), o,, is 130 meters, a standard deviation that is 10% of the 
value of the uncorrupted ocean bottom data. In both cases, the overall trend of the 
discrete bathymetric data is fit by the orthogonal function expansion. Note in Figures 
16(a) and 16(b) that the 7th-order mean-squared error and the 0th-order mean-squared 
error are of the same order of magnitude. Although a 7th-order fit is produced, the Oth- 


order, flat ocean bottom, predominates. 
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CASE: 8X10 FLAT BOTTOM WITH 1% DEVIATION IN DEPTHS 

SIGMAYB: 13.0M 7TH-ORDER MMSE SURFACE FIT 

MSFO: 1.128%10' MSE2: 1.068%10' MSE4: 1.028%10° MSE6: 1.02510" 
MSEi: 1.100«10' MSE3: 1.058%10° MSES: 1.028%10° MSE7: 1.009%10° 


Figure 16(a). Flat ocean bottom with evenly spaced data corrupted by a zero-mean 
random gaussian noise added to the depth of each data point. 
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CASE: 8X10 FLAT BOTTOM WITH 10% DEVIATION IN DEPTHS 

SIGMAYB: 130.0M 7TH-ORDER MMSE SURFACE FIT 

MSO: 1.128%10° MSE2: 1.068x10° MSE4: 1.028«10° MSE6: 1.025%10° 
MSE1: 1.100x10° MSE3: 1.068x10° MSES: 1.028x10° MSE7: 1.009%10° 


Figure 16(b). Flat ocean bottom with evenly spaced data corrupted by a zero-mean 
random gaussian noise added to the depth of each data point. 
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Figure 17 shows the shoaling ocean bottom presented in Figure 8 with Equation 


(2.24) modified to 
y,(x,z) = 1350- 3? +o,,N(@O), (2.30) 


where o, =20 meters. Again, the orthogonal function expansion fits the trend of the 


discrete bathymetric data. The values of the mean-squared errors in the legend of Figure 
17 indicate that a /st-order fit dominates, demonstrating that the orthogonal function 
expansion smoothes through noise in the discrete bathymetric data. 

More complicated ocean bottom geometries are presented in Figures 18 through 
23. The models presented in each case continue to fit the trend of the discrete 
bathymetric data and provide a smooth, continuous model of the ocean bathymetry from 
discrete data. However, the poor performance illustrated in Figure 23 demonstrates 
serious shortcomings with the two-dimensional orthogonal function expansion based on 
the current choice of generating functions. | 

Figure 18 shows a shoaling ocean bottom with a non-constant slope that increases 
as x and z increase. The exact equation used to generate the 5x 6 rectangular matrix of 


discrete data for this test case is given by 
y, (x,2) = 2000 — 3(x +z)’. (2.31) 


Equation (2.31) is of the same form as the orthogonal function expansion with N, =2. 
The mean-squared errors shown in the legend of Figure 18 show that the resulting fit is 
predominantly a 2nd-order fit. The small value of the 2nd-order mean-squared error 
shows close agreement between the discrete bathymetric data and the 2nd-order estimated 
fit. Again, the algorithm computes higher order terms, but these terms have little effect 


on the overall fit. 
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CASE: 6X8 UP-SLOPE IN Z DIRECTION ONLY, WITH WEIGHTING 

SIGMAYB: 20.0M 7TH-ORDER MMSE SURFACE FIT 

MSEO: 5.955x10° MSE2: 1.629x10* MSE4: 1.445x10° MSEG: 1.34410" 
MSE1: 1.683%10' MSE3: 1.476%10° MSES: 1.440%10° MSE7: 1.33110 


Figure 17. Planar ocean bottom sloping upward in the z direction only with no variation 
in the x direction, using W, = 0 and W, = 1, with evenly spaced data corrupted by a zero- 
mean random gaussian noise added to the depth of each data point. 
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CASE: SX6 QUADRATIC UPWARD SLOPING OCEAN BOTTOM 

SIGMAYB: 0.0M 7TH-ORDER MMSE SURFACE FIT 

MSEO: 1.022™10° MSE2: 1.803107 MSE4: 1.789%10° MSEG: 1.754%10™ 
MSE1: 7.730x10' MSE3: 1.789x107 MSES: 1.754=10° MSE7: 1.720x10° 


Figure 18. Shoaling ocean bottom in both the x and z directions 
- with evenly spaced data. 
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Figures 19 and 20 show a rolling or oscillating ocean bottom. The exact equation 


used to generate the data for Figure 19 is given by 
y,(x,z) = 1500 — 180(x +z) + 30(x +2)? -13(x +2)’, (2.32) 
and 


y,(x,z) = 1500 —180(x +z) + 30(x +z)? -13(x+z)> +o, N(0,1), (2.33) 
b a4.) 


is used to generate the data for Figure 20. In Equation (2.33), a standard deviation of 


o,, =40 meters is used to add zero-mean random gaussian noise for the test case shown 


in Figure 20. Equations (2.32) and (2.33) are of the same form as the orthogonal function 
expansion with N, =3. The mean-squared errors shown in the legend of F igures 19 and 
20 show that the resulting fits are predominantly 3rd-order. The small value of the 3rd- _ 
order mean-squared error shows close agreement between the bathymetric data and the 
3rd-order estimated fit. Figure 20 demonstrates the ability of the orthogonal function 
expansion to fit the trend of discrete bathymetric data in the presence of noise, even for | 
complex ocean bottom geometries. 

Figures 21(a) and 21(b) present results from the same test case, but shown 
graphically from two different points of view. The test case presented shows a flat ocean 
bottom with a gaussian shaped mound centered at the origin. The exact equation used to 


generate the 8 x 10 rectangular matrix of data is given by 


x? 2? 


y,(X,Z) =¢, —¢, @ 207 @ 207. (2.34) 


V2x0 
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CASE: SX6 OSCILLATING OCEAN BOTTOM 

SIGMAYB: 0.OM  7TH-ORDER MMSE SURFACE FIT 

MSEO: 2.831x10° MSE2: 2.655%10° MSE4: 1.104«10° MSEG: 1.06210" 
MSF1: 2.659x10° MSE3: 1.169%10° MSES: 1.076%10° MSE7: 1.061™10° 


Figure 19. Rolling or oscillating ocean bottom with evenly spaced data. 
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CASE: SX6 OSCILLATING OCEAN BOTTOM 

SIGMAYB: 40.0M 7TH-ORDER MMSE SURFACE FIT 

MSFO: 2.880%10° MSE2: 2.683%10° MSE4: 4.883%10° MSE6: 4.64910" 
MSE1: 2.683%10° MSE3: 4.889%10' MSES: 4.749%10° MSE7: 3,831 %10° 


Figure 20. Rolling or oscillating ocean bottom with evenly spaced data corrupted by a 
zero-mean random gaussian noise added to the depth of each data point. 
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With the constants set to c, = 1300, c, = 10°, and o = 4, Equation (2.34) becomes 


x? +27 


y,(x,Z) = 1300-994.72e 2 . (2.35) 


Equation (2.35) is modified by adding a zero-mean random gaussian noise with o, = 50 


meters, for the test case presented in Figures 22(a) and 22(b). 

Figures 21(a) and 21(b) present a test case that does not match the form of 
Equation (2.1) for the orthogonal function expansion. Although the orthogonal function | 
expansion does not exactly fit the data, the model fits the trend of the data. The gaussian 
function given in Equation (2.34) is centered at the origin. Therefore, the orthogonal 
function expansion, which is derived from the generating functions, f,(x,z) =(x+z)", 
for n= 0,1,...,7, can follow the data as it rolls away from the origin. In Figures 22(a) and 
22(b), the addition of randomness to the data creates a more realistic test case of a 
downward sloping, complex ocean bottom. The model smoothly fits through the discrete 
bathymetric data with some data points above and some data below the fit. 

Figures 23(a) and 23(b) demonstrate the major shortcoming of the algorithm 
_ being tested. This test case presents a gaussian shaped ocean bottom with the peak of the 
gaussian function centered within the rectangular matrix of data at x =5km, z=5km. 
The exact equation used to generate the 8 x 10 rectangular matrix of discrete bathymetric 
data is given by 


(x-5)? _ (2-5)? 


e 2 @ 2a | (2.36) 





y,(%,Z) = c, —C, J2no 
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CASE: 8X10 FLAT BOTTOM WITH A GAUSSIAN MOUND 

SICMAYB: O.OM 7TH-ORDER MMSE SURFACE FIT 

MSEO: 5.440«10° MSE2: 3.245*10° MSE4: 2.428™10° MSE6: 2.131™10° 
MSE1: 1.65810" MSE3: 2.819x10° MSES: 2.147%10° MSE7: 2.12810" 


Figure 21(a). Downward sloping ocean bottom with a gaussian shape in both 
the x and z directions with evenly spaced data. 
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MSE1: 1.658%10° MSE3: 2.819%10° MSES: 2.147™10° MSE7: 2.128%10° 


MSEO: 5.440%10° MSE2: 3.245%10° MSE4: 2.428%10" MSE6: 2.13110" 


CASE: 8X10 FLAT BOTTOM WITH A GAUSSIAN MOUND 
SIGMAYB: 0.0M 7TH-ORDER MMSE SURFACE FIT 


Figure 21(b). Downward sloping ocean bottom with a gaussian shape in both 
| ar 


e x and z directions with evenly spaced data. 
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CASE: 8X10 FLAT BOTTOM WITH A GAUSSIAN MOUND 

SIGMAYB: 5SO.0M 7ZTH-ORDER MMSE SURFACE FIT 

MSEOQ: 5.640%10° MSE2: 4.226%10° MSE4: 3.660%10° MSE6: 3,395x10° 
MSE1: 1.602«10° MSE3: 3.806%10° MSES: 3.39910" MSE7: 3.35310" 


Figure 22(a). Downward sloping ocean bottom with a gaussian shape in both the x and z 
directions with evenly spaced data corrupted by a zero-mean random gaussian noise 
added to the depth of each data point. 
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MSEO: 5.640%10° MSE2: 4.22610° MSE4: 3.660%10° MSE6: 3.39510" 
MSE1: 1.602™10° MSE3: 3.806%10° MSES: 3.399%10° MSE7: 3.35310 


CASE: 8X10 FLAT BOTTOM WITH A GAUSSIAN MOUND 
SIGMAYB: 50.0M 7TH-ORDER MMSE SURFACE FIT 


Figure 22(b). Downward sloping ocean bottom with a gaussian shape in both the x and 
z directions with evenly spaced data corrupted by a zero-mean random gaussian noise 
added to the depth of each data point. 
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CASE: 8X10 GAUSSIAN SHAPED OCEAN BOTTOM 

SIGMAYB: 0.0M 7TH-ORDER MMSE SURFACE FIT 

MSEO: 5.350%10° MSE2: 3.442x10° MSE4: 3.254™10° MSE6: 3.243x10° 
MSEl: 4.359x10° MSE3: 3.255x10° MSES: 3.25110" MSE7: 3.243x10° 


Figure 23(a). Gaussian shaped ocean bottom with the peak of the gaussian 
| function at x = 5km, z= 5km with evenly spaced data. 
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Figure 23(b). Gaussian 
function 
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SIGMAYB: O.O0M 7TH-ORDER MMSE SURFACE FIT | 
MSEOQ: 5.350x10° MSE2: 3.442™10° MSE4: 3.25410° MSEG: 3.243%10° 
MSE1: 4.359x10° MSE3: 3.255™10° MSES: 3.251™%10° MSE7: 3.243%10° 


CASE: 8X10 GAUSSIAN SHAPED OCEAN BOTTOM 


shaped ocean bottom with the peak of the gaussian 


at x. = 5km, z= 5km with evenly spaced data. 
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With the constants set to c, = 1300 5 C, = 10°, and o = 4, Equation (2.36) becomes 


_[(x-5)? +(2-5)7] 


y, (x,z) = 1300 — 994.72e oe 1d (2.37) 


The fit produced by the orthogonal function expansion in Figures 23(a) and 23(b) 
smoothes through the peak of the gaussian shape and provides a poor fit to the discrete 
data. This shortcoming may be due to two reasons. First, if the orthogonal function 


expansion was computed to a higher order using N, > 7, then the surface fit would more 


closely match the curvature of the gaussian shape of the discrete data. The maximum 


value of N, is a designer’s choice. The chosen value is a trade-off between exactly 


fitting rough or rapidly varying ocean bottoms versus smoothing through some of the 
roughness or rapid changes. A higher value of N, may be more appropriate if this 
technique is used to model rough or rapidly varying ocean bottoms. Second, the set of 
generating functions used to derive the equations for the two-dimensional orthogonal 
function expansion was selected due to the success of the set of generating functions used 
in the one-dimensional orthogonal function expansion presented in [Ref. 2]. The set of 
generating functions used in this thesis research may not have been the best choice. 

As can be seen in Figures 16 through 23, the surface fits have a constant value for 


the estimated ocean bottom depth across lines of contour given by 


x+z=C, (2.38) 


where C is a positive constant. The set of generating functions, f,(x,z) =(x+z)”, for 
n= 0,1,...,7, result in the same estimated depth being computed for values of x and z that 
satisfy Equation (2.38). Therefore, the surface fits have no curvature along these 
contours. This effect can be seen in both Figures 21(a) and 23(a). The estimated surface 


fits fail to roll off in two comers, the corner given by the maximum cross-range and 
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minimum down-range location, and the corner given by the minimum cross-range and 
maximum down-range location. This shortcoming is not as noticeable in Figure 21(a) 
since the discrete data has much less curvature across these contours. In Figure 23(a), 
this shortcoming results in an unacceptable fit of the discrete bathymetric data. 
Correcting this shortcoming is offered as an area of future research and possible 
alternatives will be discussed in Chapter V. 

3. Summary 

The method of orthogonal function expansion to fit two-dimensional, discrete 
ocean bathymetric data performs well, but more research is needed to fit all arbitrarily 
shaped ocean bottom geometries. The orthogonal function expansion has many positive 
features. It needs only discrete bathymetric data that may be unevenly spaced. It is 
unaffected by noise in the discrete bathymetric data, and it smoothes through noisy or 
rapidly varying data to provide a smooth, continuous model of the ocean bottom. Chapter 
III will further demonstrate the utility of this method of ocean bottom modeling by using 
this technique in the RRA Algorithm to produce accurate three-dimensional ray traces 


with rays reflecting off of arbitrarily shaped two-dimensional ocean bottoms. 
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Il. THE BATHYMETRIC MODEL APPLIED TO 
RAY TRACING IN THE RRA ALGORITHM 


A. REFLECTION ANGLES OFF THE TWO-DIMENSIONAL OCEAN BOTTOM 
MODEL | 


In a ray acoustics ray tracing program, an important feature of a model of the 
ocean bathymetry is that it facilitates accurate ray reflection calculations. The ocean 
bottom model must provide accurate first order derivatives in order to calculate the 
propagation vector of a sound ray reflecting off the ocean bottom. Additionally, in order 
for the ray tracing program to allow for three-dimensional ray propagation, a two- 
dimensional ocean bottom model is needed which can be used to calculate, or accurately 
estimate, first order partial derivatives of the ocean bottom depth with respect to both 
cross-range and down-range at the point of incidence of a propagating ray. Given these 
partial derivatives, the reflected ray can be calculated using vector analysis. | 

The partial derivatives of the ocean bottom depth with respect to cross-range and 
down-range can be accurately estimated with a central differencing algorithm that was 
originally presented in [Ref. 7], and was tested for accuracy in estimating reflection 
angles off one-dimensional ocean bottom models in [Ref. 4]. The reflection angles 
calculated in [Ref. 4] were found to be accurate, and had an average percentage error of 
less than 0.1%. The equation for estimating the partial derivative of the ocean bottom 


depth with respect to down-range is given by 


_ By, (%,Z) _ —y,(x,2+2h) +8y,(x,z +h) —8y,(x,z—-h) +y, (x,z— 2h) 


ee Aael 
oh | 12h on 


f; 


and the equation for estimating the partial derivative of the ocean bottom depth with 
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respect to cross-range is given by 


_ Oy, (x52) _ ~y,(x +2h,z)+8y,(x +h,z) -8y,(x —h,z) + y, (x - 2h,z) 


>. 43.2 
x 12h ve) 


ti 


where /; is some arbitrarily small number, with value nominally between 0.01 and 1. For 
the computer simulation test cases performed in this thesis research, a value of h = 0.1 
was used. 

Although the partial derivatives with respect to cross-range and down-range could 
be found by algebraically differentiating the equations for the orthogonal function 
expansion, the central differencing algorithm has been maintained as the method for 
estimating the partial derivatives. This technique of estimating the partial derivatives 
given any mathematical model of the ocean bottom has allowed the mathematical 
modeling technique to be changed and improved during research without having to 
change the method of calculating the partial derivatives. 

The problem geometry is illustrated in Figure 24. The angles of propagation 
shown in Figure 24(a) are taken from [Ref.1], where the direction of the unit vector along 
a ray path, 7, is determined by the angles ¢ and #. The angle # is measured from the 
positive Y axis to the unit vector, and the angle g is measured from the positive X axis to 


the projection of the unit vector into the XZ plane. The unit vector 7 can be expressed as 


AR=uUrtvyp+wz , (3.3) 

where 
u=sin Bcos¢, (3.4) 
v=cosf, (3.5) 
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and 
w=sin fsing, (3.6) 


are the dimensionless direction cosines with respect to the _X, Y, and Z axes. 
Figures 24(b) and 24(c) illustrate a sound ray reflecting off the ocean bottom. 
Figure 24(b) shows the tangent plane to the ocean bottom at the point of incidence of a 


sound ray. The incident sound ray is denoted by the unit vector 7,, and the reflected 
sound ray is denoted by the unit vector 7,. The tangent plane is defined by the unit 
normal vector to the ocean bottom, denoted 7,. Figure 24(c) presents a two-dimensional 


view of the three unit vectors shown in Figure 24(b). 
With the partial derivatives of the ocean bottom depth with respect to cross-range 


and down-range at the point of incidence of a sound ray available, the normal vector to 


the ocean bottom, denoted NV , , can be calculated as follows [Ref. 8]: 


N,=-f,8+9-f2. (3.7) 


Jeg (3.8) 





where the magnitude of the normal vector to the ocean bottom is 


W,| = Vf titf. (3.9) 


The unit vector along a sound ray at the point of incidence on the ocean bottom is denoted 
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(b) (c) 


Figure 24. (a) The unit vector along a ray path, (b) Three-dimensional view of a sound 
ray reflected off the ocean bottom in terms of the unit normal to the ocean bottom at the 
point of incidence, (c) Two-dimensional view of the vectors shown in (b). 
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by 
A, =uXtv~+w,Z. | (3.10) 


The unit vector of the reflected ray at the point of incidence, 7,, and the unit 
vector n,, have the same component perpendicular to the unit vector normal to the ocean 


bottom, 7,, but the parallel components are opposite. To obtain 7, from i,, this parallel 
component, denoted N ,, must be subtracted twice [Ref.6]. See Figures 24(b) and 24(c). 


The equation for 7, is given by 
fi, =A, -2N, =A, -2(A, 7, )A,, — B.11) 


where the inner product, (fi, -7,), is a scalar denoted NJNB. Substituting Equations (3.8) 
and (3.10) into Equation (3.11) yields 





r 


I Saag (3.12) 


A, =uUxXtv~twZz=uxtvyrwZ- NIN WW, Wy, 


where the dimensionless direction cosines associated with 7, are 





u, =U, +2 i ; (3.13) 
b 
age = 

vy, =v,—-2NINB Wy, : (3.14) 


and 
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Ww, =W, +2NINBZz, (3.15) 


Mel 


In terms of the angles # and ¢ shown in Figure 24(a), angles of reflection are 


B, =cos'(v,), (3.16) 
and 
b, = wn . (3.17) 
Uu, 


B. COMPUTER SIMULATION RESULTS 


The reflection angle algorithm presented in Section A was incorporated into the 
RRA Algorithm FORTRAN ray tracing program. Seven of the ocean bottom models 
presented in Section B of Chapter II were selected to test the accuracy of computing first 
order partial derivatives and reflection angles for propagating sound rays that reflect off 
the ocean bottom. The computer generated estimates of the partial derivatives and the 
reflection angles are compared with theoretical values. An isospeed ocean was assumed 
in each case in order to facilitate interpretation of the results. 
The seven ocean bottom models employed were: 
(1) Flat ocean bottom with evenly spaced data, as presented in Figure 2 
(2) Planar ocean bottom sloping downward in the z direction only, with no 
variation in the x direction, using W, = 0 and W, =1, with evenly spaced 
data, as presented in Figure 11 
(3) Planar ocean bottom sloping upward in both the x and z directions with evenly 
spaced data, as presented in Figure 4 
(4) Planar ocean bottom sloping upward in the z direction only, with no variation 


in the x direction, using W, = 0 and W. = 1, with evenly spaced data, as 
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presented in Figure 8 | 
(5) Planar ocean bottom sloping upward in the z direction only, with no variation 
in the x direction, using W, = 0 and W, = 1, with unevenly spaced data, as 
presented in Figure 9 
(6) Rolling bottom with evenly spaced data, as presented in Figure 19 
(7) Flat ocean bottom with evenly spaced data corrupted by a zero-mean random 
gaussian noise added to each data point, as presented in Figure 16(a) | 
In Figures 25 through 31, one ray is launched that propagates in the cross-range 
and down-range directions, with several reflections off the ocean bottom models. The ray 
is launched from an initial (x,y,z) coordinate of (0, y,,0), and initial launch angles £, 
and ¢,. Each plot shows the propagating ray’s path, a projection of the ray path onto the 
ocean bottom, and a projection of the ray path onto the ocean surface. Each plot contains 
a legend with the following information: 
(1) ‘CASE:’, a brief description of each test case 
(2) ‘Y0:’, the launch depth, y,, in meters 
(3) ‘DLTS:’, the are length step size, As , in meters, used in the Recursive Ray 
Acoustics (RRA) algorithm 
(4) ‘BETAO:’, the initial launch angle, 7, , in degrees 
(5) ‘PHIO:’, the initial launch angle, ¢, , in degrees 
For each of the seven test cases, a table of data accompanies each figure. The table 
compares computer simulation results and theoretical results for the bottom reflections. 
The theoretical values of the partial derivatives of the ocean bottom depth with respect to 


HY, (x,2Z) aga OY, (x,2) 
Cz 


cross-range and down-range at the point of incidence, —~-—— , are found 


ox 
by taking the partial derivatives of the exact equations that were used to generate the 


discrete data for each test case. The theoretical values for the angles of reflection off the 


, and the angles of the incident 


ocean bottom are computed using aan cee 
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ray 8, and ¢,. These values are applied to Equations (3.7) to (3.17) to compute 8, and 


¢,. The estimated partial derivatives of the ocean bottom depth with respect to cross- 


, are computed 


range and down-range at the point of incidence, B(x) and One 


in the RRA algorithm using the central differencing algorithm of Equations (3.1) and 


(3.2). The estimated angles of reflection for a ray reflected off the ocean bottom 


computed in the RRA algorithm are listed as B. and d,. 

Figure 25 shows a ray propagating and reflecting off the flat ocean bottom 
depicted in Figure 2. The exact fit of the ocean bottom model to the discrete bathymetric 
data results in accurate partial derivatives of the ocean bottom depth with respect to cross- 
range and down-range. For all six of the ocean bottom reflections shown in Figure 25, 


the central differencing algorithm of Equations (3.1) and (3.2) returns values of 
Seles) = SoG = —0,663172 x 10° in the computer simulation. These values on a 
good double-precision computer approximation of the actual slope of zero, that is, 
aed = a =(Q. Table 1 shows that the computer simulation eeUMe for the 


angles of propagation (in degrees) of the reflected ray are extremely accurate. 


[BOORCET aaere[ 
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Table 1. Angles of reflection off a flat ocean bottom with evenly spaced bathymetric data. 
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Figure 25. A sound ray propagating and reflecting off a flat ocean 


bottom with evenly spaced bathymetric data. 
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Figure 26 shows a ray propagating and reflecting off a planar ocean bottom 
sloping downward in the z direction only, with no variation in the x direction, as depicted 


in Figure 11. The value of the actual slope of the discrete bathymetric data with respect 


to down-range is aah = x = (0.05. The first order partial derivative calculated by 


the RRA algorithm for each bottom reflection is exactly equal to the slope of the discrete 


data, that is, Oe! = 0.05. The value of the actual slope of the discrete bathymetric 


a =(Q. The first order partial derivative 


data with respect to cross-range is 
calculated by the RRA algorithm for each bottom reflection is exactly equal to the slope 

. ; D, (x,Z) = | ae _ 
_of the discrete data, that 1s, a aa 0, due to the weighting term W, =0. Table 2 


shows that the double-precision computer simulation results for the angles of propagation 
of the reflected ray are extremely accurate. Note that although the ray is launched with an 


initial launch angle ¢, equal to 40 degrees, the ray reflections off the downward sloping 


plane cause the overall ray path to bend such that ¢, approaches 90 degrees. 
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Table 2. Angles of reflection off a planar ocean bottom sloping downward 
in the z direction only, with no variation in the x direction. 


62 





Ie 


J1G0°OF :OIHd s5900°0f :OuULIE 
WO"S =S170 WO*OOOT :0A 
NO NOILOBNIC Z NI 3dO1S-NMOG 8X8 :3SY9 


— 
(WI i 


SOILSNOOW Ave SATSHNOsY 


a planar ocean bottom sloping 


with no variation in the x direction, with evenly spaced 


Figure 26. A sound ray propagating and reflecting off 


downward in the z direction, 


bathymetric data. 
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Figure 27 shows a ray propagating and reflecting off a planar ocean bottom 
sloping upward in both the x and z directions, as presented in Figure 4. The exact fit of 
the ocean bottom model to the discrete bathymetric data results in accurate partial 
derivatives of the ocean bottom depth with respect to cross-range and down-range. For 


all six of the ocean bottom reflections shown in Figure 27, the partial derivatives 


produced by the RRA algorithm are One) = Be) = ~0,00625, which are exactly 


equal to the partial derivatives of the discrete bathymetric data. Table 3 shows that the 
angles of propagation of the reflected rays produced by the RRA algorithm are equal to 


their theoretical values. 
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Table 3. Angles of reflection off a planar ocean bottom sloping upward 
in both the x and z directions. _ 












| Figure 28 shows a ray propagating and reflecting off a shoaling ocean bottom as 
depicted in Figure 8. The first order partial derivative of the ocean bottom depth with 
respect to down-range calculated by the RRA algorithm for each bottom reflection in 


double precision accuracy is Sot) = —0.166667 , which is equal to the exact slope of 


the discrete data, that is, oS = -<. Table 4 shows that the double-precision 
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eflecting off a planar ocean bottom sloping 
with evenly spaced bathymetric data. 


Figure 27. A sound ray propagating and r 
upward in both the x and z directions, 
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RECURSIVE RAY ACOUSTICS 
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CASE: 6X8 UP-SLOPE IN Z DIRECTION ONLY, WITH WEIGHTING 
YO: 3S50.0M NUTS: 1.0” 
BETAO: S2.SDEG PHIO: 70.0DEG 
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opagating and reflecting off a planar 


ocean bottom sloping 
with evenly spaced 


Figure 28. A sound ray pr 


with no variation in the x direction, 


_ upward in the z direction, 


bathymetric data. 
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computer simulation results for the angles of propagation of the reflected ray are 


extremely accurate. 
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Table 4. Angles of reflection off a planar ocean bottom sloping upward in the z direction 
only, with no variation in the x direction, with evenly spaced bathymetric data. 








Figure 29 shows a ray propagating and reflecting off a shoaling ocean bottom 
with the discrete bathymetric data entered at unevenly spaced intervals as depicted in 
Figure 9. Table 5 shows that small errors are introduced by rounding off the input 
bathymetric data to five significant digits. The ocean bottom surface fit produced for the 
test case shown in Figures 9 and 29 is 7th-order and not /st-order due to the round off 
error introduced in the discrete bathymetric data. Overall, the effect on ray propagation ts 
small. After propagating the sound ray for a path length of over 4 kilometers, the 
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Table 5. Angles of reflection off a planar ocean bottom sloping upward in the z direction 
only, with no variation in the x direction, with unevenly spaced bathymetric data. 








67 


RECURSIVE RAY ACOUSTICS 


1250 625 | 0.2625 


is 


CASE: 6X8 UP-SLOPE IN Z ONLY, WITH UNEVENLY SPACED DATA 
YO: 350.0 DLT: 107 | 


BETAQ: 32.5DEG PHIO: 70O.0DEG 


1.0500 


cean bottom sloping 
with unevenly spaced 


Figure 29. A sound ray propagating and reflecting off a planar o 


upward in the z direction, with no variation in the x direction, 


bathymetric data. 


68 





coordinate (x,y,z) of the point of incidence of the fifth bounce in Figure 29 has been 
- displaced 2.449 meters. as compared to the coordinate (x,y,z) of the point of incidence of 
the fifth bounce of the sound ray propagated in Figure 28. Table 6 shows the theoretical 
partial derivatives of the ocean bottom depth with respect to down-range and cross-range, 
and the estimated partial derivatives computed using the central differencing algorithm at 


each bottom reflection location. 


sear | oT 


Table 6. Partial derivative of ocean bottom depth with respect to down range for 
reflections off a planar ocean bottom sloping upward in the z direction only, with no 
variation in the x direction, with unevenly spaced bathymetric data. 








Figure 30 shows a ray propagating and reflecting off a rolling or oscillating 
bottom with evenly spaced bathymetric data, as presented in Figure 19. Although the 
ocean bottom surface fit is 7th-order, the ocean bottom model closely approximates the 
discrete bathymetric data generated by Equation (2.32), which is 3rd-order. Table 7 
shows that the reflection angles produced by the RRA algorithm are accurate estimates of 
the theoretical values. Table 8 compares the estimated partial derivatives produced by the 
RRA algorithm and the partial derivatives obtained by differentiating Equation (2.32). 

The test cases shown in Figures 28 and 30 are good examples of the insight into 
ray propagation provided by using a two-dimensional ocean bottom model in the RRA 
ray tracing program. Reflections that turn the ray causing a cross-range displacement can 


be easily computed and plotted. It is important to note that the modeling of ray 
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Figure 30. A sound ray propagating and reflecting off a rolling or oscillating ocean 


bottom with evenly spaced bathymetric data. 
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Table 7. Angles of reflection off a rolling ocean bottom with unevenly spaced 
bathymetric data. 
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Table 8. Partial derivatives of ocean bottom depth with respect to cross range and down 
range for reflections off a rolling ocean bottom with unevenly spaced bathymetric data. 













propagation is intended to be an exact representation of ray propagation in the ocean. 
However, the implementation of mathematical algorithms in a computer simulation 
requires approximations to be made. The effort is to model the ocean as accurately as 
possible in order to simulate the ocean. The algorithms used in a computer model © 
designed for real-time ray tracing must be accurate and as fast and as simple as possible 
to provide fast calculations. The two-dimensional ocean bottom model produced by the 


orthogonal function expansion provides these features. 


71 





Finally, Figure 31 demonstrates the ability of the orthogonal function expansion to 
smooth through ocean bottom roughness or noisy input bathymetric data. Figure 31 
shows a ray propagating and reflecting off a flat ocean bottom with a zero-mean random 
gaussian noise added to the discrete bathymetric data as depicted in Figure 16. With a 
standard deviation of the noise at 1% of the ocean bottom depth, the ocean bottom model 
provides a smooth fit through the noise. The data in Table 9 is provided only to illustrate 
the smoothing effect of the orthogonal function expansion, and not to validate the 
reflection angle calculations. Table 9 presents data that compares the theoretical angles 
of reflection of the propagating ray reflecting off the ocean bottom for the case presented 
in Figure 25 (no noise) versus the estimated angles of reflection produced by the smooth 


fit in Figure 31 (with noise). Table 9 reports the discrepancies as a ‘% difference’ since 


BOUNCE 
8 25 a 31 | Difference | Fig. 25 | Fig.31 | Difference 


4 
5 
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Table 9. Angles of reflection off a flat ocean bottom with rough or noisy bathymetric 
data (Fig.25) versus a flat ocean bottom with no irregularities (Fig.31). 





the difference is. due to ocean bottom roughness or noisy data. The average ‘% 
difference’ between f, and £, is 0.27%, and the average ‘% difference’ between @, 


and g. is 0.20%. The error produced by successive reflections is cumulative. After 


propagating the ray over 14 kilometers, the largest difference introduced by the ocean 
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bottom roughness or noise is a 0.52% difference for the angle of reflection, B. , off the 


rough ocean bottom of Figure 31 as compared to the flat ocean bottom of Figure 25. The 


roughness or noise is effectively smoothed by the orthogonal function expansion. 
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IV. TWO-DIMENSIONAL ORTHOGONAL FUNCTION 
EXPANSIONS OF SOUND-SPEED DATA 


A. ORTHOGONAL FUNCTION EXPANSION AND THE BASIS FUNCTIONS 


Sound speed data in the ocean is typically measured as a function of depth at a 
particular geographic location. In terms of the Cartesian coordinate system presented in 
Chapters II and III, sound speed is typically measured at a fixed cross-range and down- 
range location (x,z), and the sound speed is a function of depth (y) only. This type of 
sound-speed measurement is called a sound-speed profile (SSP). The sound-speed 
measurements used in this research as test cases are discrete, and the discrete 
measurements may be taken at unevenly spaced depths. The one-dimensional orthogonal 
function expansion discussed in [Ref. 2] was demonstrated to be an effective method to 
model sound speed when it is a function of one variable, that is, depth. For a ray 
propagated down-range using the RRA algorithm, the sound-speed model developed with 
the one-dimensional orthogonal function expansion does not vary as a function of down- 
range. The one-dimensional orthogonal function expansion models sound speed as a 
function of depth only. Multiple sound-speed profiles can be used to model the speed of 
sound as a function of depth in areas (or zones) of the ocean. As the sound ray passes 
from one zone to another zone, the orthogonal function expansion of the next set of 
discrete sound-speed profile data is computed. The technique of zones is presented in 
[Ref. 2], and is illustrated in Figure 32(a). A technique is needed that will use multiple 
sound-speed profiles taken at different down-range locations to compute the speed of 
sound as a function of both depth and down-range. The two-dimensional orthogonal 
function expansion technique can be used to compute sound speed as a function of both 
depth and down-range given a rectangular matrix of discrete sound-speed data. The 
rectangular matrix of data is made up of multiple discrete sound-speed profiles taken at 
different down-range locations. The rectangular matrix of discrete sound-speed data is 


illustrated in Figure 32(b). 
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Zone | , Zone 2 | Zone 3 
sound speed | sound speed | sound speed 
inthis zone : inthiszone <: inthis zone 
isbased on : isbasedon =‘: is based on 
the SSP : the SSP : the SSP 
takenat 71 =: takenat*2 : takenat 73 

(a) 





(b) 


Figure 32. (a) Speed of sound as a function of depth (y) only, employed in 
areas (zones) of the ocean. (b) Rectangular matrix of discrete sound-speed 
data from multiple sound-speed profiles. 
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1. Definition of the Two-Dimensional Orthogonal Function Expansion 

The technique being employed in this chapter to model sound speed as a function 
of depth and down-range is the same mathematical technique that was presented in 
Chapter II to model ocean bathymetry as a function of cross-range and down-range. The 
two-dimensional orthogonal function expansion used to fit a smooth surface to a 


rectangular matrix of sound-speed data as a function of depth and down-range is given by 


é(y,2) = 16,9, (2), (4.1) 


n=0 


where the coefficient 


c, = (6, (Vy 2y)» Pe, (Yi; Zi \, 


M.. 


z 


= Cn 92) Po, VysZy) » n=0,1,---,N, ’ (4.2) 


represents the inner product of the measured discrete sound-speed data c,,(y,,2,), for 


i=1,2,...,.M 


. and =f =1,2,...,M_, with the set of orthonormal basis functions 
P, (¥.z), n=01,...,N. , evaluated at the discrete depth ( y,) and down-range(z, ) values 


corresponding to the discrete sound-speed data. The order of the expansion is represented 


by N. . The discrete data must be entered in a matrix, where M,,_ is the total number of 
depth rows and M. is the total number of down-range columns. The matrix is, in 


general, rectangular, not square, and the data is, in general, unevenly spaced in both the y 
and z directions. | | 
2. The Basis Functions from the Generating Functions 


As in Chapter II, the set of orthonormal basis functions ?,. (y,Z), n= 0,1,...,M.; 


were determined by using the Gram-Schmidt procedure [Ref. 3] on the set of generating 
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functions f,(y,z)=(y+z)", n=0,1,...,N,. The generating function chosen is of the 
same form as in Chapter II. The general equation of the Gram-Schmidt 


Orthogonalization Procedure is given by 


(F,(9,2).0.,(02)) 
(0, (¥3Z),P., (y.z)) 
(F,(1:2)s@.,, 2) 


ap Le, (2), n=012,.... (4.3) 
(9, (y,Z), PQ. (y,z)} 


Po, (y;2) = Ff, (V2) < P., (.2)=%: , 


The Gram-Schmidt procedure to find the basis functions, 9, , for n=0,1,...,N,, is the 


same as the procedure detailed in Section A of Chapter II, and it will not be repeated here. 
As a means of comparison to the derivation shown in Chapter II, the final equations for 
the first three basis functions are listed below: 
The Oth-order basis function is given by 
P., (y,2) = I : (4.4) 
The /st-order basis function is given by 
P,.(¥.z)=yV+z—-Wy . (4.5) 
The 2nd-order basis function is given by 


2. (y.z)=(y +z) - Wy -Waly+z-Wo) . (4.6) 


The use of the two-dimensional orthogonal function expansion to model sound 


speed as a function of depth and down-range gives similar results as the model of 
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bathymetric data presented in Chapter I]. The data may be unevenly spaced. As in 


Chapter II, where N, =7 is the maximum order fit, the maximum order sound-speed fit 


is N.=7. Noise in the input discrete sound-speed data does not affect the overall fit. 


_ Also, the mean-squared error (MSE) is computed for NV, = 0,1,...,7. The value of N. 


that produces the lowest mean-squared error is used to reconstruct the sound-speed data. 
B. COMPUTER SIMULATION RESULTS OF SOUND-SPEED FITS 


The two-dimensional orthogonal function expansion technique to model sound- 
speed data was implemented in the RRA Algorithm FORTRAN ray tracing program. 
The results of various test cases are presented graphically to demonstrate both the 
successes and the shortcomings of this method of modeling discrete sound-speed data. 
The discrete sound-speed data is made up from three sound-speed profiles (sound-speed 
versus depth) taken at three different down-range locations. The graphs, presented in 
Figures 33 to 41, show the discrete sound-speed data points as large black dots (e). The 
estimated fit to the discrete data computed by the orthogonal function expansion is 
presented as a mesh of lines to represent the two-dimensional, continuous surface. 

The technique of using a two-dimensional orthogonal function expansion to 
provide a continuous surface fit to discrete data is the same technique that was tested in 


Chapter Il. The general conclusions and shortcomings presented in Chapter I apply in 


this chapter. Even with the shortcomings discussed in Chapter II, the testing of this 


technique to fit sound-speed data will validate the FORTRAN subroutines that compute 
the two-dimensional orthogonal function expansion of sound-speed data in the RRA 
algorithm computer program. The test cases will also provide specific conclusions and 
shortcomings of this technique when it is applied to fit discrete sound-speed data. 

Each plot contains a legend with the following information: 

(1) ‘CASE’, a description of the test case 

(2) ‘SIGMAC:’, the value of the standard deviation, o ,, of a randomly 
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generated zero-mean gaussian noise added to each sound-speed data point 
(3) the order of the fit, N. , based on the minimum mean-squared error 
criterion 
(4) ‘MSEO:’ to ‘MSE7:’, the mean-squared error for NV, = 0,1,2,...,7 
Figures 33 and 34 demonstrate the fit of the most basic test case, a constant speed 
of sound. The exact equation used to generate the 6x3 rectangular matrix of data for 


Figures 33 and 34 is given by 
c(y,z) = 1500. (4.7) 


In both test cases, the fit is an exact Oth-order fit with the orthogonal function expansion 
returning a value of 1500 m/s for any combination of depth and down-range. Figure 34 
shows that the data may be entered at unevenly spaced coordinates without affecting the 
fit. | 

Figures 35, 36, and 37 demonstrate exact /st-order fits. The exact equation used to 
generate the 6 x 3 rectangular matrices of data for Figures 35, 36, and 37 is given by 


c(y,z) = 1490+ 0.02y, (4.8) 


which represents sound speed that is linearly increasing as a function of depth. Figure 35 
reveals the same shortcoming that was found with the /st-order bathymetric fits presented 
in Chapter II. The corner of the surface fit at the maximum depth and minimum down- 


range location (y = 1000m, z = 0.0 km) is below the discrete sound-speed data points; the 


comer of the surface fit at the minimum depth and maximum down-range location 


(y = 0.0m, z= 1.0km) is above the discrete sound-speed data points. The fitted surface 


approximates a /st-order planar surface that slopes upward equally in both the y and z 


directions. Therefore, the fitted plane is below half of the discrete sound-speed data 
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CASE: CONSTANT SOUND SPEED WITH EVENLY SPACED DATA 
SIGMAC: O.0M OTH-ORDER MMSE SURFACE FIT 

MSEQ: 4.098x10° MSE2: 4.09810. MSE4: 4.09810° MSEG: 4.098%x10™ 
MSEL: 4.098™10° MSE3: 4.098%10° MSES: 4.098%10° MSE7: 4.098x10™ 


Figure 33. Constant sound-speed with evenly spaced data. 
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CASE: CONSTANT SOUND SPEED WITH UNEVENLY SPACED DATA 

SIGMAC: 0.0M OTH-ORDER MMSE SURFACE FIT 

MSEQ: 4.098x10°% MSE2: 4.098x10° MSE4: 4.098%10° MSE6: 4.098%10™ 
MSE1: 4.098%10°% MSE3: 4.098x10° MSES: 4.098%10° MSE7: 4.09810" 


Figure 34. Constant sound-speed with unevenly spaced data. 
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CASE: SS LINEARLY INCREASING, EVENLY SPACED DATA NO WT 
SIGMAC: O.0M  7TH-ORDER MMSE SURFACE FIT | 
MSEO: §.400%107 MSE2: 4.941%10° MSE4: 4.827%10° MSES6: 4.825*10° 
MSE1: 4.941%107 MSE3: 4.827107 MSES: 4,825»10° MSE7: 4.66310" 


Figure 35. Sound-speed that is linearly increasing as a function of depth, 
with no variation in the down-range direction, modeled without 
weighting x and z separately, with evenly spaced data. 
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CASE: SS LINEARLY INCREASING, EVENLY SPACED DATA WITH WT 

SIGMAC: 0.0M 1ST-ORDER MMSE SURFACE FIT 

MSEQ: 8.400%10° MSE2: 5.215x10° MSE4: 5.215=10° MSE6: 5.21510" 
MSE1: 5.215x10° MSE3: 5.215x10° MSES: 5.215=10° MSE7: 5.21510" 


Figure 36. Sound-speed that is linearly increasing as a function of depth, 
with no variation in the down-range direction, using W, =1 and W, =0, 


with evenly spaced data. 
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SOUND SPEED 





CASE: SS LINEARLY INCREASING, WITH UNEVENLY SPACED DATA 
SIGMAC: 0.0M §7TH-ORDER MMSE SURFACE FIT 

MSFO: 7.895%10° MSF2: 1.321«10% MSE4: 1.311%10° MSEG: 1.275x10° 
MSE1: 1.370102 MSE3: 1.317107 MSES: 1.290107 MSE7: 1.157x10° 


Figure 37. Sound-speed that is linearly increasing as a function of depth, 
with no variation in the down-range direction, using W, =1 and W, = 0, 


with unevenly spaced data. 
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points in the lower half of Figure 35, and above half of the discrete sound-speed data 
points in the upper half of Figure 35. The best fit is limited in its accuracy due to the fact 
that y and z are summed in every term of the orthogonal function expansion. As in 
Chapter II, an independent weighting with respect to y and an independent weighting with 
respect to z need to be added to the orthogonal function expansion to correct this 
shortcoming. With the addition of weighting factors, the orthogonal function expansion 


becomes 


N. 
é(y,z) = 06,9, (W,yW,2) (4.9) 


n=Q 


where W, and W, are weighting factors to introduce independence in y and z. 


Figure 36 presents the same test case presented in Figure 35, except that the 
weighting factors, W, and W,, were added as shown in Equation (4.9) with W, =1 and 
W.=0. The result is an exact /st-order fit. Figure 37 presents the same test case with 
the discrete sound-speed data entered at unevenly spaced locations. The fit determined by 
the orthogonal function expansion is 7th-order due to using only five significant digits for 
the discrete sound-speed data. The overall fit is /st-order; the higher order terms of the 
orthogonal function expansion only affect the lower significant digits of the estimated 


value of sound speed. 


Figures 38, 39, and 40 present a test case with discrete sound-speed data 


generated by the exact equation given by 


c(y,z) = 1525-0.05y, (4.10) 


which represents sound-speed that is linearly decreasing as a function of depth. In Figure 
38, no weighting is applied, and the orthogonal function expansion produces a poor fit. 
In Figures 39 and 40, the weighting factors, W, and W,, were added as shown in 
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CASE: SS LINEARLY DECREASING, EVENLY SPACED DATA NO WT 
SIGMAC: 0.0M § 7TH-ORDER MMSE SURFACE FIT 

MSEQ: 5.250x10° MSE2: 3.088«10° MSE4: 3.017%10° MSE6: 3.016™10° 
MSE1: 3.088%10° MSE3: 3.017x10° MSES: 3.016™10° MSE7: 2.914x10° 


Figure 38. Sound-speed that is linearly decreasing as a function of depth, 
with no variation in the down-range direction, modeled without 
weighting x and z separately, with evenly spaced data. 
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SOUND SPEED 





CASE: SS LINEARLY DECREASING, EVENLY SPACED DATA WITH WT 
SIGMAC: O.0M 1ST-ORDER MMSE SURFACE FIT 

MSEO: 5.250x10° MSE2: 7.078%10° MSE4: 7.078%10° MSE6: 7.07810" 
MSE1: 7.078™x10° MSE3: 7.078%10° MSES: 7.078™%10° MSE7: 7.078™10™ 


Figure 39. Sound-speed that is linearly decreasing as a function of depth, 
with no variation in the down-range direction, using W, = 1 and W_ = 0, 


with evenly spaced data. 
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CASE: SS LINEARLY DECREASING, WITH UNEVENLY SPACED DATA 

SIGMAC: 0.0M | ZTH-ORDER MMSE SURFACE FIT | 
MSEO: 5.120*10° MSE2: 7.723x10° MSE4: §.346x10° MSE6: 6.023%10™ 
MSE1: 8.196x107 MSE3: 7.66810 MSES: 6.024x10° MSE7: 6.013x10™ 


Figure 40. Sound-speed that is linearly decreasing as a function of depth, 
with no variation in the down-range direction, using W, = 1 and W_ =0, 


with unevenly spaced data. 
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Equation (4.9) with W,=1 and W,=0. An exact /st-order fit results in Figure 39. 


Figure 40 demonstrates that the discrete sound-speed data may be unevenly spaced 
without affecting the overall fit. 

| Figure 41 demonstrates a more complicated test case. The discrete sound-speed 
data used in this test case represents three different sound-speed profiles taken at three 
different down-range locations. The discrete data for the constant sound-speed profile at 
z= 0.0km is generated by Equation (4.7). The discrete data for the linearly decreasing 
sound-speed profile at z = 0.5km is generated by Equation (4.10). The discrete data for 


the parabolic sound-speed profile at z = 1.0km is generated by the exact equation given 


by 
c(y,z) = 1490 + 0.0001(y — 500)’. (4.11) 


The three sound-speed profiles form a 6x3 rectangular matrix of discrete sound-speed 
data. The result from the orthogonal function expansion is a smooth surface fit that 
attempts to transition from the constant sound-speed profile to the linearly decreasing 
sound-speed profile to the parabolic sound-speed profile. The smooth transition is 
obtained. However, at the down-range locations where the sound-speed profiles were 
taken, the characteristics of the speed of sound as a function of depth have been distorted. 
The distortion can be seen more clearly by examining a two-dimensional plot of the fitted 
speed of sound versus depth taken at the fixed down-range locations of the three sound- 
speed profiles. 

Figures 42, 43, and 44 show two-dimensional plots of the speed of sound versus 
depth for the test case presented in Figure 41. The plots show the discrete sound-speed 
data points as large black dots (¢ ). The estimated fit to the discrete data computed by the 
two-dimensional orthogonal function expansion is presented as a continuous curve. 


Significant distortion resulted. 
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CASE: VARYING SCUND SPEED, CONSTANT-LINEAR-PARABOLIC 

SIGMAC: O.0M 7TH-ORDER MMSE SURFACE FIT 

MSEQ: 2.358«10° MSE2: 2.157%10° MSE4: 9.409x107 MSE6: 8.813%10° 
MSE1: 2.182«10° MSE3: 1.109%10° MSES: 8.929%107 MSE7: 8.45710" 


Figure 41. Sound-speed that varies in the down-range direction, with the discrete - 
sound-speed data constant as a function of depth at z = 0.0 km, linearly decreasing as a 
function of depth at z = 0.5 km, and parabolic as a function of depth at z = L.0km, 
modeled without weighting x and z separately, with evenly spaced data. 
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0.0 km for the test case shown in Figure 41, 


ed versus depth at z = 
with the discrete sound-speed data constant as a function of depth at z 
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Figure 42. Sound- 


0.0 km, modeled 


without weighting x and z separately, with evenly spaced data. 
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Figure 43. Sound-speed versus depth at z 


0.5 km for the test case shown in Figure 41, 


with the discrete sound-speed data linearly decreasing as a function of depth at 
0.5 km, modeled without weighting x and z separately, with evenly spaced data. 
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CASE: VARYING SOUND SPEED, CONSTANT-LINEAR-PARABOLIC SSP AT Z = 1.000 KM 

ZONE: 1 SIGMAC: O.OM/S 7TH-ORDER MMSE SURFACE FIT 
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Figure 44. Sound-speed versus depth at z = 1.0km for the test case shown in Figure 41, 


with the discrete sound- 
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The distortion of the fitted sound-speed as compared to the discrete sound-speed 
profiles is due to the smoothing effect of the orthogonal function expansion. If the 


orthogonal function expansion was computed to a higher order using N. > 7, then the 


sound-speed surface fit would more closely match the changes in the discrete sound- 
speed profiles. The distortion is also a result of the set of generating functions used. As 
discussed in Chapter II, the surface fits have a constant value for the estimated sound 


speed across lines of contour given by 
y+z=C, | (4.12) 


where C is a positive constant. Therefore, the surface fits have no curvature along these 
contours. The amount of distortion for different sound-speed profiles will depend on how 
well the rectangular matrix of discrete sound-speed data matches the form of the 
generating functions. For a given set of data, the fit that is obtained may be acceptable, 
but the technique is not versatile enough to provide acceptable fits in all cases. 

As in Chapter II, the method of orthogonal function expansion to fit two- 
dimensional, discrete data provides a starting point for the use of this technique to 
compute a continuous model of the speed of sound as a function of depth and down- 
range. The orthogonal function expansion needs more variability in the maximum order 
fit that is selected, N,. The technique could be improved by using a different set of 
generating functions that will not return constant values for the speed of sound across 
lines of contour given by Equation (4.12). For two-dimensional sound-speed fits, the use 
of this technique has an additional shortcoming. For the sound-speed profiles used as 
input data, their basic shape as a function of depth may be distorted. A better and simpler 
technique may be to use a one-dimensional orthogonal function expansion to fit each 
sound-speed profile separately, and then average or interpolate between the nearest one- 


dimensional sound-speed fits. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This thesis had two primary goals. First, to derive and test the two-dimensional 
orthogonal function expansions of ocean bathymetric and sound-speed data. This 
technique provides a continuous mathematical model, from unevenly spaced discrete 
input data, of bathymetry as a function of cross-range and down-range, and the speed of 
sound as a function of depth and down-range. The second goal was to implement the 
expansions in the Recursive Ray Acoustics (RRA) Algorithm FORTRAN program to 
evaluate the successes and shortcomings of this technique to model the ocean bottom and 
the speed of sound. In order to implement the technique as the ocean bottom model in the 
RRA Algorithm, a method of computing accurate angles of reflection also had to be 
developed. The orthogonal function expansion has many positive features. It takes 
discrete input data that may be unevenly spaced, and provides a fit to the discrete data that 
smoothes through noise or rapid variations in the data. The result is a continuous two- 
dimensional model of the discrete data. The major shortcoming is that the technique 1s 
limited by the choice of the set of generating functions used. Also, the maximum order 
fit that is implemented has been limited to seventh order. The choices for the set of 
generating functions and the maximum order fit were based upon previous research with 
the one-dimensional orthogonal function expansion [Ref. 2]. 

The effort to provide a more sophisticated one-dimensional ocean bottom model 
for the RRA Algorithm, other than a flat ocean bottom, began in 1991 with the use of 
Spatial Fourier Series [Ref. 4]. The technique met with limited success. The one- 
dimensional orthogonal function expansion was developed and presented in 1994 [Ref. 
2]. The one-dimensional orthogonal function expansion was successful, and the 
technique was used to model both bathymetric and sound-speed data in the RRA 
Algorithm. The two-dimensional orthogonal function expansion may yet prove to be a 
viable tool in modeling both ocean bathymetric data and sound-speed data, but its 


immediate utility, based on the analysis conducted, 1s limited. 
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It is recommended that future research apply higher-order fits and test generating 
functions that will correct the shortcomings produced by the set of generating functions 
used in this analysis. A possible alternative technique could be a two-dimensional 
orthogonal function expansion that is based on the generating functions of the one- 
dimensional case [Ref. 9]. The form of the orthogonal function expansion used to fit a 
rectangular matrix of discrete ocean bathymetric data as a function of cross-range and 


down-range would then be given by 


N, Ns 


By (2.2) = DD Pam Ps, Ys, (2), (5.1) 
n=0 m=0 
where 9, (x), 2=0,1,...,.V,, and y, (z), m=0,l,...,N,, are the orthogonal functions 


obtained by applying the Gram-Schmidt Orthogonalization Procedure to the generating 


functions given by 

f,(x)=x", n=0,1...,N,, (5.2) 
to find g, (x), and 

g,(zZ)=z", m=0,l,...,N,, | (5.3) 


to find y, (z). The coefficients 5,,,, 2=0,1,...,N, and m=0,l,...,N,, would be 
obtained by computing the inner product of the measured discrete ocean bathymetric data 
with the orthonormal basis functionsg, (x)y, (z). This technique would provide a 
greater degree of freedom to fit the surfaces with more independence in x and z. For a 


-seventh-order fit, N,=7, this technique would produce 64 coefficients 
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(b.95p.15-+-207.6677) as opposed to the eight coefficients developed in the technique 


presented in Chapter II. 

The ray trace plots presented in Chapter III provide insight into the complexity of 
ray propagation in the ocean. The study of sound propagation involving multiple ocean 
bottom reflections requires a good two-dimensional ocean bottom model. With the 
continued interest of the United States Navy in understanding the acoustic environment in 


the littorals, more advanced study in this technique of mathematical modeling is 


warranted. 
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